arXiv:astro-ph/0512063vl 2 Dec 2005 


Mon. Not. R. Astron. Soc. 000. H UTUl f2002'! Printed 5 February 2008 (MN LMjTjX style file v2.2) 

Imaging and integral field spectroscopy of shocked H 2 
around G25.65+1.05 

Stephen P. Todd 1,2 * and Suzanne K. Ramsay Howat 1 f 

1 UK Astronomy Technology Centre, Royal Observatory, Blackford Hill, EH9 3HJ, UK. 

2 Institute of Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, EH9 3HJ, UK. 

Accepted . Received ; in original form 


ABSTRACT 

Near-infrared imaging of the emission from molecular hydrogen is a powerful 
method for discovering outflows in star-forming regions. We present new near-infrared 
images, long slit and integral field spectroscopy of the ultra-compact H II region 
G25.65+1.05. These new observations reveal shocked H 2 emission associated with a 
bipolar outflow from a young high mass star at the centre of the source. The physical 
parameters of the outflow are discussed and compared with outflows from lower mass 
stars. 
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1 INTRODUCTION 

A number of recent surveys have invigorated the study 


of high mass star 
herent samples of c 
iects (HMYSOs) (e.e 

formation by providing data on co- 
andidate high mass voung stellar ob- 

.. Palla et alJ 1991, Molinari et al.1 1996, 

Molinari et all 1998, 

Ridge & Moore 200ll Sridharan et alJ 

120021). Until these : 

surveys became available the relative 


scarcity of HMYSOs, the large distances to the closest ex¬ 
amples (a few kpc) and the high extinction to these ob¬ 
jects had limited their study to a few, well-known examples. 
Thus, there are many questions still remaining as to the na¬ 
ture of high-mass star formation, and the similarities and 
differences in high and low mass stellar evolution. 

Outflows are associated with the formation of 
stars of all masses and are the subject of _a_ number 
of the HMYSO survey s (IShepherd^fcjChurchwehl Il996bl 
IShepherd fc ChurchwelA ll996al . lzhang^UrHT20^T) . While 
the role of outflows in the formation of low-mass stars is 
com parativ ely well studied (see for example the review by 
iRicher et al.l 1200(1) . far less is known about the properties 
and role of outflows associated with the formation of high- 
mass stars (M greater than ~ 8 Mg). A radio survey of 
molecular line emission from high-mass star forming regions 
showed that high-velocity molecular gas (CO ) is associated 
with around 90% of these regions dShepherd fc Churchwelll 
ll996a ). Two sources (G25.65+1.05 and G240.31+0.0T) were 
mapped at higher spatial resolution and found to have bipo- 
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lar outfl ows. In a follow-up survey, fshepherd fc Churchweli 
Jl99f>bh mapped a further ten of the best candidates for hav¬ 
ing bipolar molecular outflows and confirmed th e presence 
of such an outflow in 5 of those sources. IZhang et ahltiflOll 
observed 69 candidate high mass protostars in the CO J=2-l 
transition and argue that as many as 90% of the sources may 
have outflows. iBeuther et a 0 12002a), mapping at a higher 
spatial resolution, found evidence of bipolar outflows in 21 of 
their 26 sources suggesting that bipolar outflows may be as¬ 
sociated with most young high mass stars. The first report of 
molecular hydrogen line e missi o n as sociated with high mass 
star formation was from iLee et all ( 2001 1, following up a 
search for the near-infrared count e rparts of ultra-co mpact 
H II regions bv lWalsh et alJ dl999ll . ILee et, all ll200 ill found 
evidence for outflows in IRAS 15278-5620 and IRAS 16076- 
5134. 

These surveys confirm the importance of outflows for 
HMYSO evolution. Further work is required to determine 
whether the physics of the HMYSO outflows is the same 
as for low mass YSOs. The observed outflows appear to 
have much lower collimation factor - between 1 and 1.8 - 
than those seen from low-mass stars which often have a col¬ 
limation factor of around 10 (see, for examp le, the inter¬ 
ferometric observations of IRicher et al.1 (,' 20001 1. This would 
be hard to explain if the outflows are formed by the same 
jet entrainment model as that believ ed to describe t he out - 
flows from low-mass stars. However. iBeuther et a D <l2002alt 
argued that the observed degree of collimation could be sig¬ 
nificantly reduced by the low spatial resolution of the maps, 
which could be consistent with well collimated high-mass 
flows. Interferometric observations have shown collimation 
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factors as high as 10 dBeuther et al.112002b ') in flows from 
high-mass stars and have revealed that some of these appar¬ 
ently uncollimated flows can be resolved into several well col¬ 
limated flows from separate young stars. Detection of mul¬ 
tiple flows in a region in which high-mass stars are forming 
would not be sur prising. High-mass star s are known to form 
in dense clusters llGarav fc Lizandl999i) . so high spatial res¬ 
olution observations are essential for identifying the source 
of an individual outflow. The presence of a collimated out¬ 
flow would imply the presence of a stable accretion disk and 
hence strengthen the view that high mass stars are formed 
by steady accretion, in a similar way to low mass stars, 
rather than by merging of intermediate-mass protostars in 
the centre of dense clusters. 

We have embarked upon a near-infrared survey of 
HMYSOs, which will be presented in full in Varricatt et 
al. (2005, in preparation). The survey is designed to reveal 
outflows by mapping at sub-arcsecond spatial resolution in 
the v=l-0 S(l) line of molecular hydrogen. NIR observa¬ 
tions of H 2 have frequently been used to study outflows 
from low mass YSOs. In this paper we present near-IR imag¬ 
ing, integral-held and long-slit spectroscopic observations of 
new outflows in the region of the ultra-compact H 11 region 
G25.65+1.05 . 

G25.65+1.05 (also IRAS 18316-0602 or RAFGL7009S) 
is an irregular, com p act radio source, first identified at 
3.6cm byjKurt^^hl 1 W9$). It is located at a distance of 
3.2kpc liMolina^e^d] 1996 1. The radio peak is coincident 
with an unresolved infrared source, identihe d as a young 
B1V star with a large K-band excess llZavagno et al.| |2002l) 
and i s also closely associated with methanoljjMolinane^y 
ll996|) and ammonia maser emission ifWaMm^d J •ioo.l 
ISzvmczak et al ] 200(1). Submillimetr e continuum observa- 
tions at 350p m ||Hunt er et alJ l2000l) . 450pm and 850pm 
llWalsh et alJ12003) " are all peaked at the position of the 
ra dio and maser sources. Observation of the CS (2-1) line 
bv lBronfman et’ahl (1996;) shows an excellent match of the 
observed ra dial velocity from the masers (40.8-42.4 kms -1 , 
IWalsh et al.11200 3) and the line emission (41.4 kms -1 ) in¬ 
dicating a strong link between the dense gas, the m aser 
sources and the massive star. IZavagno et, alJ 120021 pro¬ 
pose that one explanation for the observed A-band ex¬ 
cess from the central source may be the presence of a 
disk. ISO spectroscopy of this source shows a rich spec¬ 
trum of ice features including absorption features attributed 
to H 2 0, CH 3 OH, C0 2 , 13 C0 2 ,CO, PC S, HCOOH HCO Q-, 
CH 3 HCO, CH 4 , NH 3 and Silicate (see lGibb et all2004l and 
references therein). Laboratory spectra fitted to the ISO ob¬ 
servations suggests that the ice features arise in dense ma¬ 
terial with temperatures in the range 10K-100K. 

2 OBSERVATIONS 
2.1 Imaging 

Imaging observations were made on 29^June 2002 using the 
facility infrared imager, UFTI dRoche et al.h2QQ 3). at the 
UK Infrared Telescope (UKIRT). UFTI is a 1 to 2.5 pm 
camera using a 1024 x 1024 HgCdTe array with a plate scale 
of 0.091 arcsec/pixel. A single 9-point jitter in a 3 x 3 pattern 
with offsets of 10 arcsec was observed with 60 sec exposures 




Figure 1. Part of the UFTI mosaic images (a) K98 broad-band 
filter; (b) 1 — 0 S(l) narrow-band filter. All coordinates are J2000. 


using the A98 broad-band filter, giving 9 min on source. 
The same 9-point jitter pattern was repeated three times 
using exposures of 100 s with the 2.122 pm 1 — 0 S(l) H 2 
narrow-band filter giving a total of 45 min on source. The 
individual frames from each filter were flat-fielded and mo- 
saiced together automatically by the ORAC-DR pipeline. 

Sources emitting in the H 2 line will be detected in both 
the broad-band and narrow-band images but will appear 
brighter relative to the continuum sources in the narrow- 
band images. The two mosaics are shown in Figure 1. The 
narrow-band image with 0.6 arcsec seeing was smoothed to 
match the 0.7 arcsec seeing of the broad-band image. The 
broad-band image was scaled in intensity using the flux from 
the field stars and subtracted from the narrow-band image 
to leave the H 2 emission as shown in Figure 2. The H 11 
region itself and some stars around it appear negative in 
this subtracted image due to differential extinction to the 
core compared with the field stars used for flux scaling. The 
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Figure 2. A continuum subtracted 1 — 0 S(l) image created by subtracting Figure la from Figure lb. 
The contours are drawn from a version of the image smoothed with a 0.5 arcsec gaussian and are 
logarithmically spaced with each contour a factor of two higher in flux than the previous one from 
an arbitrary lowest contour. All coordinates are J2000. 


increased reddenening makes it relatively brighter in the 
broad-band image. A number of point sources are visible 
in addition to the extended emission. These are ghost im¬ 
ages of bright stars produced by the narrow-band filter and 
residuals left where stars have been imperfectly subtracted. 
The complete mosaic covers a 1.8 x 1.8 arcmin field centred 
on the H II region. No H 2 emission was detected outside the 
region shown here. The brightest H 2 emission was detected 
to the south-east of the H 11 region, taking the form of a re¬ 
gion of faint, diffuse emission containing a number of bright, 
compact sources (A-D). To the north-west of the H 11 re¬ 
gion there is a straight, narrow line of faint emission with 
a bright source (E) at one end. There is also faint, diffuse 
emission to the south and north-east of the H 11 region. 

2.2 Long-slit spectroscopy 

Long slit spectra of sources C and D were obtained on 30 
June 200 2 at UKIRT using Cooled Grating Spectrometer 4 
llMountain et a HH). The 40 1/mm grating was used, giv¬ 
ing a wavelength coverage of 1.9 to 2.5 pm and the sources 
were observed with a single slit position. A one pixel (0.6 arc- 
sec) wide slit was used and the detector was stepped over 2 
pixels in half-pixel increments to give a fully sampled spec¬ 
trum with a spectral resolution of 600. Observations were 


made through thick, patchy cirrus causing large variations 
in the detected flux. Individual exposures were weighted by 
the square of the signal to noise ratio of the v = 1 — 0 S(l) 
line before they were added together to maximise the signal 
to noise ratio of the combined data. The spectrum was ra- 
tioed by a standard star with its blackbody shape corrected, 
to remove atmospheric telluric features. The final spectra of 
the two sources are shown in Figure 3. 

2.3 Integral field spectroscopy 

Sources A and B were observed on 24 October 2002 as part of 
the UIST commissioning observations. UIST is a new facility 
class near-IR (1-5 pm) imager and sp ectrometer at UKI RT 
which uses a 1024 x 1024 InSb array llR.amsav Howat et alJ 
l2000il . It includes a deployable cryogenic image-slicing IFU 
with 14 slices, each 0.24 arcsec wide. The pixel scale along 
the slices is 0.12 arcsec/pixel resulting in a field of view 
of 3.36 arcsec x 5.52 arcsec, which can be rotated to any 
angle on the sky. The HK grism was used, giving a spectral 
coverage of 1.4-2.5 pm with a spectral resolution of 800- 
1000 . 

The target was acquired using UIST in A-band imag¬ 
ing mode. The IFU field of view was rotated to a position 
angle of 90°, making the long axis of the field run E-W. Two 
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a) 



b) 



Figure 3. CGS4 long-slit spectra of a) source C and b) source 
D. The flux is in arbitrary units. 



Figure 4. The spectral line image of sources A and B at the 
wavelength of the 1—0 S(l) line, extracted from the IFU datacube. 
The offsets are arbitary and from the edge of the IFU field, which 
was centred on a = 18 h 34 m 23.5 s , 5 = — 6 o 53'00 / '. At a distance 
of 3kpc, larcsec corresponds to 0.015pc. 


adjustments to the pointing were made during the observa¬ 
tions (1 arcsec E, then 1 arcsec S). A total of three positions 
were observed for 8mins each giving a maximum integra¬ 
tion of 24mins on-source for the overlap region. Sky frames 
were obtained by offsetting to a separate sky position. The 
standard star BS7260 was observed to provide correction for 
telluric atmospheric features and flux calibration. 

2.3.1 Combining the data into a single datacube 

The observations at each of the three positions described 
above were formed into an (x, y, A) data-cube using the 
ORAC-DR pipeline employed for all on-line data reduction 
at UKIRT. Each frame was divided by a flat-field frame. 
Sky-subtracted frames were wavelength calibrated using a 
spectrum of an argon arc lamp and the resulting frames 
formed into a datacube. All of the spectra in the datacube 
were divided by the standard-star spectrum, corrected for 
the black-body shape of the stellar continuum, to remove 
variations in atmospheric transmission. An approximate flux 
calibration was obtained using the magnitude of the star 
(K=4.48). The IFU data reduction pipeline is described in 
more detail in lTodd et alJ (l2002h . 

The known 1 arcsec telescope offsets were used to reg¬ 
ister the three data-cubes in the two spatial dimensions. 
These were mosaiced into a single data-cube using make- 
MOS, which is included in the Starlink CCDPACK package. 
An image extracted from this datacube at the wavelength 
of the v = 1 — 0 S(l) H 2 line (2.122 pm) is shown in Figure 4. 
A spectrum of the sum of the brightest part of source A is 
shown in Figure 5. All detected emission lines are from H 2 . 

A spectrum was formed by summing over all spatial 
pixels of the data-cube. A Gaussian profile was fitted to the 


bright 1 — 0 S(l) line in this spectrum and was found to 
be centred at (2.12209 ± 0.00003) pm with FWHM (2.28 ± 
0.05) x 1CU 3 pm. The other lines are all much weaker, so 
to obtain the best possible measurement of the flux in each 
line a Gaussian profile was scaled to fit the line keeping the 
FWHM of the line fixed at 0.00228 pm and the offset of the 
line from the 1 —0 S(l) line fixed by the known wavelength of 
the line. In total, 12 lines from the v=l and v=2 transitions 
were detected. The measured fluxes are reported in Table 1. 

3 RESULTS AND DISCUSSION 
3.1 Morphology 

The K-band image of G25.65+1.05 (Figure la) is dom¬ 
inated by the nebulous continuum emission surrounding 
the central, unresolved, source at a = 18 h 34 m 20.9 s ,8 = 
—6°02'42.3" seen by Zavagno et al. (2002). The morphol¬ 
ogy of the H 2 emission in Figure 4 is strongly suggestive of 
a moderately well-collimated outflow, or outflows, centred 
on G25.65+1.05. Sources A, B, C and D are connected by 
diffuse emission and appear to form one lobe of an outflow 
with a flow axis through G25.65+1.05 to source E. The out¬ 
flow has a position angle on the sky of 130° East of North. 
This is consistent with the direction of the highly energetic 
bipolar outfl ow fou nd to be centred on or close to the H 11 
region bv IShepherd fc Churchwefil lll996bFl . Sources ABCD 
are spatially co-incident with the blueshifted lobe seen in 
the CO outflow and E with the redshifted lobe. The rela¬ 
tive faintness of E compared with ABCD may then be in¬ 
terpreted as being due to E being embedded further in the 
molecular cloud. The projected length of the outflow, taken 
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Figure 5. The spectrum obtained by summing over the brightest region of source A. Fluxes measured from 
this spectrum are given in Table 1. 


Table 1. Measure fluxes and physical parameters for the lines measured. 


Transition 

A j / jam 

E j /K 

Aj /10“ 7 s 

9j 

Flux /W m 2 ; 

arcsec 2 

1 

- 0 S(7) 

1.7480 

12817 

2.98 

57 

(2.4 ±0.2) 

X 

10“ 16 

1 

- 0 S(3) t 

1.9576 

8365 

4.21 

33 




1 

- 0 S(2) 

2.0338 

7584 

3.98 

9 

(4.4 ±0.2) 

X 

IQ - 16 

2 

- 1 S(3) 

2.0735 

13890 

5.77 

33 

(1.7 ±0.2) 

X 

IQ - 16 

1 

— 0 S(l) 

2.1218 

6956 

3.47 

21 

(11.4 ±0.2) 

X 

10~ 16 

1 

- 0 S(0) 

2.2235 

6471 

2.53 

5 

(2.7 ±0.2) 

X 

10 -16 

2 

- 1 S(1) 

2.2477 

12550 

4.98 

21 

(1.7 ±0.2) 

X 

IQ - 16 

1 

— 0 Q(l) t 

2.4066 

6149 

4.29 

9 




1 

- 0 Q(2) t 

2.4134 

6471 

3.03 

5 




1 

- 0 Q(3) t 

2.4237 

6956 

2.78 

21 




1 

- 0 Q(4) 

2.4375 

7586 

2.65 

9 

(3.0 ±0.3) 

X 

IQ - 16 

1 

- 0 Q(5) t 

2.4548 

8365 

2.55 

33 





The wavelengths, upper energy levels, Einstein A coefficients and degeneracies of 
the H 2 lines detected in our spectra. Lines marked ' were found to be absorbed by 
unresolved atmospheric absorption features and were discarded. The fluxes given are 
measured from the spectrum in Figure 5 and are not corrected for extinction. 


to be the distance from C to E, is ~70arcsec. In the absence 
of any information, we assume an angle relative to the plane 
of the sky of 45°. This gives an estimate of the total length 
of 1.4pc and a collimation factor ~3. Interpreted as a single 
outflow, we find a lower degree of collimation than is seen 
in outflows from low-mass YSOs, but higher than the 1-1.8 
often assumed for high-mass YSOs. 

IShepherd fc Churchwelll lll996bll found two velocity 
components associated with G25.65+1.05 at 39kms~ 1 and 
at 49kms . One possible interpretation of the H 2 emission 
in this region would be to see sources A and B as bow-shocks 
in a highly collimated flow from west to east, as suggested 
by their morphology. This would imply the presence of an¬ 
other, more deeply embedded source to the south an d west of 
G25.65+1.05. The MSX point source catalogue ( Eean et al . 


{1992) shows a single source associated with G25.65+1.05 
(a = 18 h 34 m 21.2 s , <5 = —6°Q2 l 'i6.Q6" with a positional accu¬ 
racy of 4” —5”) so does not support an argument for a second 
source offset south-west of G 25.65+ 1.05 . Fu rthermore, a re¬ 
cent submillimetre survey bv lHill et alJ d200Rll finds a single 


cated at the position of the methanol maser 

Molinari et al. 

199a), consistent with an earlier result by 

Faundez et al. 


In the following Section, the near-IR spectra of sources 
ABCD are analysed and shown to support our model of the 
emission arising due to a single outflow. 
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3.2 H 2 Excitation in the G25.65+1.05 Outflow 


The near-infrared emission from H 2 can be produced ei¬ 
ther by thermal excitation in shock-fronts or fluorescent ex¬ 
citation by non-ionizing ultraviolet photons (91.2 < A < 
110.8 nm) from hot young stars. These cases may be distin¬ 
guished using the near-IR spectrum, though the detection 
of lines from a broad range of energy levels is required to 
obtain a secure result. For shock-excited H 2 , the lower en¬ 
ergy levels (v=l) are typically populated as for a gas in local 
thermal equilibrium (LTE) with a characteristic excitation 
temperature of a few thousand Kelvin, though in reality the 
temperature of the cooling gas will vary over a range from 
a few hundred to a few thousand Kelvin. The emission lines 
would also be broadened to tens of kms -1 . This is less than 
our instrumental resolution and undetectable in the spec¬ 
troscopy of G25.65+1.05. 

For radiatively excited gas, the populations follow a 
non-LTE distribution, characterised by high excitation tem¬ 
peratures (~ 10,000K). Lines from high vibrational levels 
(v>> 1 ) may be detected and the ratio of the intensity 
of the 1 — 0 S(l) line to that of the 2—1 S(l) H 2 line is 
~ 2 if the density is low llSternberg fc Dalgarndll989 b For 
gas above the critical density (nH 2 ~ 10 5 cm -3 ), collisional 
de-excitation becomes important and the level populations 
tends towards an LTE population. 

The spectra of G25.65+1.05 provide high signal/noise 
observations of lines from the v=l and v =2 vibrational lev¬ 
els. Deeper observations would be required to obtain line 
strengths for the higher vibrational levels and thereby to 
confirm the excitation. However, there are indicators that 
the H 2 is shock excited which, taken in aggregate, lead us 
to conclude that this is indeed the case. 

In the following sections, we obtain the excitation tem¬ 
peratures for the sources in G25.65+1.05. These are in the 
range 1000-3000K indicative of thermal excitation. When 
molecular hydrogen is radiatively excited, 90% of the H 2 
is dissociated by the ionizing photons (A < 91.2 nm) 
that are also present. This may result in emission in the 
2.166 pm Brackett -7 recombination line. The spectra of 
sources ABCD contain only emission lines of H 2 . We note 
that the absence of Brackett -7 emission is not a strong con¬ 
straint as some fluorescently excited photodissociation re¬ 
gions emit only H 2 lines. Finally we cite the bow-shock mor¬ 
phology seen clearly in Figure 4, perhaps the strongest in¬ 
dicator that the H 2 is excited by a shock front driven by an 
outflow rather than heated from within the knots by embed¬ 
ded point sources. 


3.2.1 The Excitation Temperature of Sources A and B 


The measured intensity, I, of a given H 2 line can be used to 
calculate the column density of the upper excitation level of 
the transition: 


Nj = 


4-7tA jl 
Ajhc 


(1) 


where Aj is the Einstein A-coefficient of the transition. The 
relative column densities of any two excitation levels can be 
expressed in terms of an excitation temperature T ex : 


Ni_ 

Nj 



f — {Ej — Ej) \ 
\ kT ex ) 


( 2 ) 


where gj is the degeneracy and Ej is the energy of the level. 
The values of Aj, Ej , Aj and gj for the lines detected in our 
spectra are shown in Table 1. 

Before we could use our measured intensities to derive 
the excitation temperature of the gas it was necessary to 
measure and compensate for extinction. In the absence of 
more information it was assumed that the extinction was 
constant across the IFU field of view. An extinction law 
of the form r(A) = Ak (A/2.2pm ) -1 ' 75 was used, giving a 
corrected intensity of 7 cor r = I/e~ T( ‘ x \ Plotting log(A/,/g_,) 
against Ej should give a straight line for low vibration levels 
and for a single characteristic excitation temperature in each 
pixel. If the value of Ak used to correct the line intensities 
from which the column density is calculated is wrong then 
the scatter of the points will increase. Our measurements are 
particularly sensitive to this because we have measured one 
A-band line (1 — 0 S(7) at 1.748 pm) which comes from an 
upper level close in energy to the upper level of two A-band 
lines (2 — 1 S(l) at 2.248 pm and 2—1 S(3) at 2.074 pm). 

The value of Ak which minimised the scatter of the 
points was measured for each spatial pixel in the two bright¬ 
est columns in the rebinned data-cube. It was found that 
the intensities of the 1 — 0 S(3) line and all the Q-branch 
lines other than the 1 — 0 Q(4) line were not consistent with 
any non-negative values of Ak- These lines were therefore 
assumed to be partly absorbed by the very narrow atmo¬ 
spheric lines which dominate the edges of the A-band win¬ 
dow which, being spectrally unresolved, are not removed 
by dividing by the standard star. Once these lines were ex¬ 
cluded the mean value of Ak over all the spatial pixels used 
was found to be Ak = 0.7 ± 0.1, where the error on the 
derived value was estimated from the scatter of the values 
measured from one spatial pixel to another. The variation 
in the measurements appeared random, with no evidence 
for any systematic variation in extinction across the source. 
There was also no evidence of curvature in the Nj /gj versus 
Ej plots, one of which is shown in Figure 6 , though the ab¬ 
sence of any points with upper energy levels between 8000 
and 12000 K would make it hard to detect small deviations 
from a straight line. A curved line in this plot would provide 
evidence for non-LTE processes, such as would be seen for 
radiatively excited gas. 

The line fluxes were measured at each spatial pixel 
in the rebinned data-cube and corrected for extinction us¬ 
ing the extinction law described above. The value of Nj /gj 
was then derived from each corrected line flux using Equa¬ 
tion (1). The excitation temperature was measured by fitting 
Equation (2) to all of the measurements made at a single 
spatial pixel using a maximum likelihood method. When a 
single intensity is used there is a degeneracy between T and 
the constant of proportionality, causing the most likely value 
of the constant of proportionality to vary exponentially with 
T. For this reason Equation (2) was reformulated as 

f =6XP ( : # +Q ) (3) 

where a is the logarithm of the constant of proportionality. 

The likelihood of a range of values of T and a was calcu¬ 
lated for each measured line flux assuming a gaussian error 
distribution on the measurements of the line intensity and 
using a uniform prior probability density function for both 
parameters: 
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Eypper /K 


Figure 6. Plotting Nj/gj against E U pper on a logarithmic scale 
gives a straight line after correcting for extinction using the ex¬ 
tinction law described in the text. The solid line is given by Equa¬ 
tion (3) using the maximum likelihood values of T and a. These 
measurements are taken from the 0.48arcsecx0.48 arcsec region 
of the source marked ’d’ Figure 4. 


Lj (T, a) oc 


exph(<y2o-) 2 ] 

a\/2n 


(4) 


In Equation (4) 5 is the difference between the value 
of Nj /gj derived from the measured intensity using Equa¬ 
tion (1) and that calculated using Equation (3) and a is the 
lcr error on the measured value of Nj /gj. When this is eval¬ 
uated for a single spectral line there is a degeneracy between 
the two parameters, as shown in Figure 7a. The likelihood 
of the parameters using all the line fluxes is calculated by 
taking the product of all the individual likelihoods: 


L(T,a) = Y[Lj(T,a) 

3 


(5) 


This reduces the degeneracy (Figure 7b). The likelihood of 
each value of T, taking into account the uncertainty in a, can 
be calculated by marginalising the likelihood, or summing 
over all values of a to produce a one dimensional likelihood 
curve (Figure 8). The 68% confidence levels (the narrowest 
possible range of values of T to include 68% of the total area 
under the likelihood curve) were measured as an equivalent 
of In error bars on a gaussian distribution. 

These measurements were carried out along three 
0.48 arcsec wide strips marked a, b and c on the image shown 
in Figure 4. The results are shown in Figure 9. The excita¬ 
tion temperature clearly increases from around 1800 K to 
(2840 ± 230) K (averaging over the most easterly pixel of all 
three slices) at what may be the shock front. In source A 
there also appears to be a decrease in excitation tempera¬ 
ture of ~ 200 K to the right (east) in the image shown in 
Figure 4. 

In Figure 4 it can be seen that both sources are brighter 
on the left (south) side than on the right. The temperatures 
measured in source A are approximately 200 K higher on 
the south edge of the IFU field than 1.5 arcsec to the north. 
We interpret this asymmetry in the north-south direction 
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Figure 7. (a) The likelihood of the parameters T and a were 
evaluated for a single spatial pixel using the intensity of a single 
spectral line, (b) When other spectral lines are introduced the 
degeneracy is partly broken. On both plots the contours are at 
the 68%, 95% and 99% confidence levels, equivalent to lcr, la and 
3(7 error bars on a gaussian distribution. 




Figure 8. The marginal distribution of T can be obtained by 
integrating the distribution shown in Figure 7b over a. The like¬ 
lihood shown here is normalised to have a maximum value of 1. 
The maximum likelihood value of T (2178 K in this case) and con¬ 
fidence intervals can be measured from this. The 68% confidence 
interval of 1973 — 2371 K is shown by the dotted lines. 


as supporting our interpretation of the H 2 emission being 
associated with a single outflow. If sources A and B are were 
produced by a flow from west to east a symmetrical bow- 
shock would be expected. However, we cannot exclude the 
possibility that the asymmetry is produced by variations in 
density in the medium with which the outflow is interacting. 











8 S. P. Todd and S. K. Ramsay Howat 


0 - 0.48 arcsec 


3000 

2500 

2000 

1500 


0 1 2 3 4 5 6 7 

Offset along slice /arcsec 


0.48 - 0.96 arcsec 



Offset along slice /arcsec 

0.96 - 1.44 arcsec 


3000 

2500 


1500 


0 1 2 3 4 5 6 7 

Offset along slice /arcsec 

Figure 9. The temperature was measured along three 0.48 arcsec 
wide strips running west-east on the left of Figure 4. The points 
and error bars on these plots show the maximum likelihood value 
and 68% confidence interval on the value of T marginalised over 
a, as described in the text. 


3.2.2 The Excitation Temperature of Sources C and D 

The spectra of sources C and D are shown in Figure 3 and 
the relative line intensities in Table 2. Although the weather 
was non-photometric, both sources were observed simulta¬ 
neously and the line ratios may be considered to be accu¬ 
rate, at the limit of the gaussian errors. The maximum like¬ 
lihood analysis was repeated for C and D, giving an exci¬ 
tation temperature of 2730 K with 68 % confidence levels of 
(2727±227)K for source C and (2256±250)K for source D, 
in each case averaged over the observed extent of the source. 


3.2.3 Shock excitation 

The properti es of four classes of shocks were summarised by 
iDavies et al.l (l2000Tl : 

(i) fast J shock (100-300 km s -1 ): hydrogen molecules are 
dissociated and reform producing a spectrum similar to that 
of UV fluorescence. The flux of Brackett -7 and 77-band Fe 11 
lines are comparable to the 1 — 0 S(l) flux. 

(ii) slow J shock-, strong H 2 lines are produced and 
molecules are not dissociated. However, in normal ioniza¬ 
tion fractions, magnetic field strengths and gas densities, 
such shocks are expected to evolve into C-type shocks. 

(iii) fast C shock: shock velocities of ~ 40 km s -1 heat 
the gas to ~ 2000 K or more (dependent on the velocity) 
producing strong H 2 lines. 

(iv) slow C shock: peak temperature 300 K or less, result¬ 
ing in very weak H 2 emission. 

The absence of Brackett -7 emission from atomic hy¬ 
drogen or Fe 11 lines and the excitation temperatures mea¬ 
sured in all the sources ABCD leads us to suggest that the 
fast C shock model may be most appropriate. A maximum 
post-shock excitation temperature of (2840 ± 230) K would 
be produc ed by a shock velocity of ( 37 ± 3) kms -1 in the 
models of lKaufman fc Neufeldl ill 9961) . This would be con- 
sis tent w ith the 29 km s " 1 alon^_the line of sight measured 
bv lShepherd fc Churchwel l (11996^ when the uncertainty on 
the angle of incidence is taken into account. The tempera¬ 
tures measured in source A (offsets along the slice of between 
4.5 and 6.5 arcsec) decrease with di stance from the shock 
front. The C-shock models of lFlower et alJ ( 1996h predict a 
steady decrease in temperature over around 0.03 pc (2 arcsec 
at a distance of 3 kpc) behind the shock front, beyond which 
point the flux falls rapidly. This is consistent with what we 
see, though the proximity of source B and the uncertain 
geometry of the outflow inhibits more detailed comparison 
with the models. A more detailed analysis, including lines 
from higher vibrational levels examined as function of dis¬ 
tance across the shock-front or high spectral resolution to 
explore the lines shapes would be required to allow us to 
draw stronger conclusions. Integral field spectroscopy is a 
technique that would lend itself readily to these observa¬ 
tions. 


4 SUMMARY 

Mapping at high spatial resolution in the near-IR. emission 
from H 2 we have directly observed outflow activity from 
one or more sources associated with G25.65+1.05. Long slit 
and integral field spectroscopy of the brightest sources show 
that they are excited by a shock. The measured excitation 
temperatures and absence of dissociation lead us to suggest 
that a C-shock with a shock velocity of (37±3) kms -1 may 
be responsible. Integral field spectra of a region containing 
two compact sources of H 2 emission showed a variation in 
excitation temperature with a maximum of 2800 K close to 
the shock front. The excitation temperature in the brighter 
of the two sources falls steadily over around 0.03 pc at which 
point the intensity dec reases significa ntly, in agreement with 
the C-shock models of lFlower et alJ (.1996). 

Both the location of the brightest H 2 emission and 
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Table 2. Line ratios from the spectra of 
sources C and D. 


Line Ratio to 1-0 S(l) 

Source C Source D 


1-0 

S(l) 


1 



1 


2-1 

S(l) 

0.15 

± 

0.02 

0.10 

± 

0.04 

1-0 

S(0) 

0.20 

± 

0.05 

0.36 

± 

0.07 

2-1 

S(3) 

0.22 

± 

0.04 

0.15 

± 

0.05 

1-0 

S(2) 

0.28 

± 

0.08 

0.52 

± 

0.08 


The ratio of line fluxes measured in the CGS4 
spectra of sources C and D relative to the 1- 
0 S(l) line flux. The spectra were observed in 
cloudy conditions, so the absolute fluxes are 
unknown. 


the derived shock speed are consistent with the excited H 2 
tracing the same outflow as t he CO emission detected by 
IShepherd fc Churchwelll | 1996 bh. 

We conclude that the most likely interpretation of these 
data is that the H 2 emission traces the full width of the 
outflow. I 11 this case, the collimation factor of the outflow is 
~ 3, higher than the 1 to 1.8 often though to be typical of 
outflows from high mass stars. 

The combination of near-infrared imaging and spec¬ 
troscopy is a powerful method of revealing outflows in re¬ 
gions of high mass star formation. However, in the case of 
G25.65+1.05, the near-infrared data alone cannot determine 
unambiguously whether there are multiple or single out¬ 
flows and driving sources. Supporting information may be 
sought in a number of ways. Higher resolution spectroscopy 
can help associated the emission with different outflows (see 
e.g. iDa'vis et, al.l (120041) on IRAS18151-1208). High resolu¬ 
tion imaging at wavelengths from 3-20microns and longer, 
such as is now possible from large ground based telescopes, 
may reveal additional deeply embedded sources. 

The combinination of high spatial resolution interfero¬ 
metric radio mappin g with inf r ared i m aging has been used 
succ essfully by e.g . jLeeet alJ ll200ll) . IWalsh et, alJ (l2002ll 
and iBeuther et alJ (12003) to examine the detailed struc¬ 
ture and morphology of crowded regions containing mul¬ 
tiple sources and outf lows. G323.7 4-0.26 is revealed by the 
observations of lWalsh et alJ (2002) to contain both fluores¬ 
cent H 2 emission from the Hu region associated with this 
source and shocked H 2 emission from an outflow. Combining 
interferomet ric observations wit h infrared images of IRAS 
19410+2336 IBeuther et al.l i2003fl found that what had ap¬ 
peared in single dish maps to be two outflows was the com¬ 
bination of at least seven and possibly as many as nine sep¬ 
arate outflows. It seems likely that similar results might be 
found in this region. 
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